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We study theoretically the optical response of the surface states of a topological insulator, espe- 
cially the generation of helicity-dependent direct current by circularly polarized light. Interestingly, 
the dominant current, due to an interband transition, is controlled by the Berry curvature of the 
surface bands. This extends the connection between photocurrents and Berry curvature beyond the 
quasiclassical approximation where it has been shown to hold. Explicit expressions are derived for 
the (111) surface of the topological insulator Bi2Se3 where we find significant helicity dependent 
photocurrents when the rotational symmetry of the surface is broken by an in-plane magnetic field 
or a strain. Moreover, the dominant current grows linearly with time until a scattering occurs, 
which provides a means for determining the scattering time. The dc spin generated on the surface 
is also dominated by a linear-in-time, Berry curvature dependent contribution. 



I. INTRODUCTION 

Topological insulators (TIs) are characterized by topo- 
logically protected surface states (SSs). In their simplest 
incarnation, these correspond to the dispersion of a sin- 
gle Dirac particle, which cannot be realized in a purely 
two dimensional band structure with time reversal invari- 
ance. This dispersion is endowed with the property of 
spin- momentum locking, i.e., for each momentum there 
is a unique spin direction of the electron. Most of the 
experimental focus on TIs so far has been towards trying 
to directly observe these exotic SSs in real or momentum 
space, in tunneling[l| and photoemissionf^, H, HI exper- 
iments, respectively, and establish their special topolog- 
ical nature. However, there has so far been a dearth of 
experiments which study the response of these materials 
to external perturbations, such as an external electro- 
magnetic field. 

In order to fill this gap, we study here the response of 
TI surfaces to circularly polarized (CP) light. Since pho- 
tons in CP light have a well-defined angular momentum, 
CP light can couple to the spin of the surface electrons. 
Then, because of the spin-momentum-locking feature of 
the SSs, this coupling can result in dc transport which is 
sensitive to the helicity (right- vs left-circular polariza- 
tion) of the incident light. This phenomenon is known as 
the circular photogalvanic effect (CPGE). In this work, 
we derive general expressions for the direct current on a 
TI surface as a result of the CPGE at normal incidence 
within a two-band model and estimate its size for the 
(111) surface of Bi2Se3, an established TI, and find it to 
be well within measurable limits. Since bulk Bi2Se3 has 
inversion symmetry and the CPGE, which is a second- 
order non-linear effect, is forbidden for inversion symmet- 
ric systems, this current can only come from the surface. 

We find, remarkably, that the dominant contribution 
to the current is controlled by the Berry curvature of the 
electron bands and grows linearly with time. In prac- 
tice this growth is cut-off by a scattering event which 
resets the current to zero. At the microscopic level, this 
part of the current involves the absorption of a photon 



to promote an electron from the valence to the conduc- 
tion band. The total current contains two other terms - 
both time-independent - one again involving an interband 
transition and the other resulting from intraband dynam- 
ics of electrons. However, for clean samples at low tem- 
peratures, the scattering or relaxation time is expected 
to be large, and these contributions will be eclipsed by 
the linear-in-time one. Hence, this experiment can also 
be used to measure the relaxation time for TI SSs. 

Historically, the Berry curvature has been associated 
with fascinating phenomena such as the anomalous Hall 
effect [3] and the integer quantum Hall effect Q and there- 
fore, it is exciting that it appears in the response here. 
Its main implication here is that is gives us a simple rule, 
in addition to the requirement of the right symmetries, 
for identifying the perturbations that can give a linear- 
in-time CPGE at normal incidence: we look for pertur- 
bations that result in a non-zero Berry curvature. Put 
another way, we can identify perturbations that have the 
right symmetries but still do not give this current because 
the Berry curvature vanishes for these perturbations. Im- 
portantly, for TI SSs, the requirement of a non-zero Berry 
curvature amounts to the simple physical condition that 
the spin-direction of the electrons have all three compo- 
nents non-zero. In other words, if the electron spin in 
the SSs is completely in-plane, the Berry curvature is 
zero and no linear-in-time CPGE is expected. The spins 
must somehow be tipped slightly out of the plane, as 
shown in Fi^QJi, in order to get such a response. Thus, 
a pure Dirac (linear) dispersion, for which the spins are 
planar, cannot give this response; deviations from linear- 
ity, such as the hexagonal warping on the (111) surface 
of Bi9Te?;(l0|. are essential for tilting the spins out of the 
plane. 

CPGE has been observed in the past in GaAsQ and 
SiGe[7] quantum wells - both systems with strong spin- 
orbit coupling. However, no connection with the Berry 
curvature was made in those cases. The connection, 
if present, may be harder to find because the descrip- 
tion there necessarily involves transitions between four 
bands - two spin-orbit split valence bands and two spin- 
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Figure 1: (Color online) (a) Schematic illustration of pref- 
erential absorption at one out of two points related by the 
reflection symmetry about the yz-plane. The short arrows 
denote the spin direction of electrons in various states. At 
low energies, the spins are completely in-plane. They acquire 
a small out-of-plane component at higher energies. The dot- 
ted lines represent incoming photons of helicity —1 (left-CP 
photons). These photons can only raise the (S z ) of an elec- 
tron, and thus are preferentially absorbed by electrons whose 
(S z ) < in the valence band. The chemical potential /j, must 
be between the initial and final states for any absorption to 
occur, (b) Constant energy contours for the surface conduc- 
tion band of Bi2Se3. Dark lines denote lower energy, (a) is 
drawn at p y = 0. (c) Geometry of the experiment. Light is 
incident normally on (111) surface of Bi2Se3. The dotted lines 
represent the mirror plane m about which the lattice has a 
reflection symmetry. The current j a 2{t) (see text) is along x. 



assuming a scattering time of lOps, (the scattering time 
in GaAs is ~ Ins over a wide range of temperatures (l3| ; 
we use a conservative estimate for Bi2Se3 here) we find 
that a current density of ~ lOOnA/mm (~ lOnA/mm) 
can be obtained due to the CPGE with a TWatt laser. 
This value can be easily measured by current experimen- 
tal techniques. Conversely, the scattering time, crucial 
for transport processes, for Bi2Se3 SSs can be determined 
by measuring the current. In comparison, circular photo- 
galvanic currents of a few nanoamperes per Watt of laser 
power have been measured in GaAs and SiGe quantum 
wells. 

A connection between the optical response of a system 
and the Berry curvature of its bands has been previously 
noted at the low frequencies, where a semiclassical mech- 
anism involving the anomalous velocity of electrons in 
a single band explains itjlij. Here, we show it for in- 
ter band transitions where no quasiclassical approxima- 
tion is applicable. Instead, we calculate the quadratic 
response function directly. A connection is still present 
which points to a deeper relation between the response 
functions and the Berry curvature. 

This paper is organized as follows. In Sec. [Ill we 
state the symmetry conditions under which a CPGE may 
occur. We present our results, both general as well as 
for Bi2Se3 in particular, in Sec. MI Al and describe the 
microscopic mechanism in Sec. 1111 B[ The calculation is 
described briefly in Sec. MI CI and in detail in Appendix 
IB"1 In Sec. IIV1 we give our results for dc spin. 



degenerate conduction bands. In contrast, TI SSs can be 
faithfully treated within a two-band model. Thus, TI SSs 
are a more convenient system theoretically compared to 
semiconductor quantum wells. In general, if a surface has 
no rotational symmetry about the surface normal, such 
a photocurrent is allowed. 

Finally, we estimate the current on the (111) surface 
of Bi2Se3 using an effective model for the SSs@, |To| . 
This model captures the deviations from linearity of the 
SS dispersion due to the threefold rotational symme- 
try of the (111) surface of Bi 2 Se 3 . These deviations 
have been observed in photoemission experiments on 
Bi2Te3 0. Similar deviations are expected for Bi2Se3Q, 
though they cannot be seen in the slightly smaller mo- 
mentum range compared to Bi2Te3 over which data is 
currently available [Ty. In order to get a direct current 
with CP light at normal incidence, rotational symmetry 
about the surface normal needs to be broken. Based on 
the requirement of non-zero Berry curvature, we propose 
to do this in two ways: 

1. by applying an in-plane magnetic field and includ- 
ing deviations from linearity of the dispersion 

2. by applying a strain. 

With a magnetic field of 10T (With a 1% strain) and 



II. SYMMETRY CONSIDERATIONS FOR THE 
CPGE 

In this section, we specify the symmetry conditions 
under which one can get a CPGE on the surface of a TI. 
But first, let us briefly review the concept of the CPGE 
in general. 

The dominant dc response of matter to an oscillating 
electric field is, in general, quadratic in the electric field. 
When the response of interest is a current, the effect is 
known as the photogalvanic effect. This current can be 
written as 

3a = 7?a^ 7 ^/3(w)f 7 (-w) (1) 

where £ a (t) = £ a (oj)e iut + £^(uj)e~ lult is the incident 
electric field, £&((*]) — £ a {~u) and ry Q /3 7 is a third rank 
tensor, which has non-zero components only for systems 
that break inversion symmetry, such as the surface of a 
crystal. 

For j a to be real, one has rj a p 7 = ?7* 7 ^- Thus, the real 
(imaginary) part of r] a ^ is symmetric (anti-symmetric) 
under interchange of /3 and 7, and therefore describes 
a current that is even (odd) under the transformation 
lu — > —u). Consequently, j a can be conveniently sepa- 
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rated according to 



£p(u)£*{u) + £* p (ui)£^(ui] 



-ik ail (£x£ ) fl 
(2) 

where S a ^ 7 is the symmetric part of r) a p 7 and A afl 
is a second-rank pseudo-tensor composed of the anti- 
symmetric part of rjaPj- For CP light. £ oc x ± iy if z 
is the propagation direction and only the second term in 
Eq. @ survives, and hence represents the CPGE. This 
effect is odd in tu. On the other hand, the first term, 
which is even in tu, represents the linear photogalvanic 
effect as it is the only contribution for linearly polarized 
light. Since the transformation uj — > — lj, or equivalently, 
£ — > £* reverses the helicity of CP light, i.e., changes 
right-CP light to left-CP light and vice versa, the CPGE 
is the helicity-dependent part of the photogalvanic effect. 

The helicity of CP light is odd (i.e., right- and left- 
CP light get interchanged) under time-reversal. It is also 
odd under mirror reflection about a plane that contains 
the incident beam, but invariant under arbitrary rota- 
tion about the direction of propagation. Let us consider 
normal incidence of CP light on a TI surface normal to 
the z axis. Let us further assume that there is a mirror 
plane which is the y-z plane (See Fig. [Th). Then the sym- 
metries above imply that the only component of direct 
current that reverses direction on switching the helicity is 
a current along the x axis. If there is also rotation sym- 
metry R z about the z-axis (such as the threefold rotation 
symmetry on the (111) surface of Bi2Se3), then no sur- 
face helicity-dependent direct photocurrent is permitted. 
One needs to break this rotation symmetry completely 
by applying, for example, and in-plane magnetic field, 
strain etc., to obtain a nonvanishing current. 



III. HELICITY-DEPENDENT DIRECT 
PHOTOCURRENT 

We now present our main results for the photocurrent 
and estimate it for Bi2Se3. After painting a simple micro- 
scopic picture for the mechanism, we give a brief outline 
of the full quantum mechanical treatment of the phe- 
nomenon. 



cr are the spin-Pauli matrices. Clearly, this can cap- 
ture a Dirac dispersion, eg. with E(p) = ±vfp and 
n(p) = vpz x p. It can also capture the SSs of Bi2Se3 in 
the vicinity of the Dirac point, which includes deviations 
beyond the Dirac limit. We also assume the Hamiltonian 
has a reflection symmetry to about y-axis, where z is the 
surface normal. Using the zero temperature quadratic 
response theory described in Sec MI CI we calculate the 
current due to the CPGE and find that 



jcPGE(t) = (jna+jal +ia2^)) X 



(4) 



where the subscripts a (na) stand for "absorptive" and 
"non-absorptive", respectively. The absorptive part of the 
response involves a zero momentum interband transition 
between a pair of levels separated by energy too. These 
terms are only non zero when there is one occupied and 
one empty level. In this part of the response, we find a 
term that is time-dependent, j a 2(t)- In particular, this 
term grows linearly with the time over which the elec- 
tromagnetic perturbation is present, which is allowed for 
a dc response. In reality, this linear growth is cut off 
by a decay process which equilibrates populations, and is 
characterized by a time constant r. In clean samples at 
sufficiently low temperatures, characterized by large r, 
this contribution is expected to dominate the response, 
and hence, is the focus of our work. The other contribu- 
tions are discussed in Appendix [Bj Conversely, because 
of the linear growth with time, one can determine the 
lifetime of the excited states by measuring the photocur- 
rent. This term is 



i«2(t) 



7re 3 ft£Qtegn(uj) 



Y,5(H\u)\-2\E p \)v x (p)F(p) 



. (5) 

where we have assumed that the chemical potential is 
in between the two energy levels ±|S P | connected by 
the optical frequency hui, and that temperature can 
be neglected compared to this energy scale. Here, 
v x(p) = dgp P is the conventional velocity and F(p) = 

iJ2p(dp* u (p)\dpy u (P)) + c - c -> wnere Kp)) is the con- 
duction band Bloch state at momentum p, is the Berry 
curvature of the conduction band at momentum p. For 
the class of Hamiltonians (|3|) that we are concerned with, 
the Berry curvature is given by (See Appendix [A]) : 



Results 



A general two-band Hamiltonian (in the absence of the 
incident light) can be written as 



£# P = £|£p|fi(p)-<r 



(3) 



upto a term proportional to the identity matrix, which 
is not important for our main result which involves only 
inter-band transitions. Here n(p) is a unit vector and 



F(P) 



dh dh 

^dp x dp y 



(6) 



which is the skyrmion density of the unit vector n in 
momentum space. Since d Pi h _L n for i — x,y, F(p) ^ 
only if all three components of n are nonvanishing. 
For linearly dispersing bands, n has only two non-zero 
components (eg. H p = p y a x - p x a y , n oc (p y , -p x , 0)). 
Hence, corrections beyond the pure Dirac dispersion are 
essential. Also, due to to, the Berry curvature satisfies 
F(p x ,p y ) = —F(—p x ,py). Since in Eq. §5§ we have the 
x-velocity multiplying the Berry curvature, which also 
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transforms the same way, a finite contribution is obtained 
on doing the momentum sum. 

We now calculate j a 2(t) for the threefold-symmetric 
(111) surface of Bi2Se3 starting from the effective 
Hamiltoniand [TjJ 

H = v F (p x a y - p y cr x ) + ^ (p+ + pl_) a z (7) 

where v F ~ 5 x 10 5 m/s[lJ] and A = 50.1eV-A 3 @|. A 
spin independent quadratic term has been dropped since 
it does not modify the answers for interband transitions, 
which only involve the energy difference between the 
bands. 

To get a non-zero jcPGE, the threefold rotational sym- 
metry must be broken, which we first propose to do by 
applying a magnetic field B in the x-direction. This field 
has no orbital effect, and can be treated by adding a 
Zeeman term —g x fJ.BBa x , where g x is the appropriate g- 
factor and \ib is the Bohr magneton, to the Hamiltonian 
d?]). To lowest order in A and B, we get 



3a2{t) 



3e 3 v F £ 2 X(g x [i B B) 2 t 
16h 2 uj 



A 



(8) 



to lowest order in A and B 2 , where A is the laser spot- 
size. For g x = 0.5[H, and assuming the experiment is 
done in a 10T field with a continuous wave laser with 
hui = O.leV which is less than the bulk band gap of 
0.35eV@], A ~ 1mm 2 , a laser power of 1W, and the 
spin relaxation time t ~ lOps, we get a current density of 
~ lOOnA/mm, which is easily measurable by current ex- 
perimental techniques. Note that the expression (JHJ| for 
ja.2 (t) contains the parameter A which measures the cou- 
pling to a z in Eq. |7]) . Since B — Bit breaks the rotation 
symmetry of the surface completely, a naive symmetry 
analysis suggests, wrongly, that deviations from linear- 
ity, measured by A, are not needed to get j a 2{t)- 

The rotation symmetry can also be broken by applying 
a strain along x, which can be modeled by adding a term 
SXpla z to H in EqEJ. This gives 



. 3eV(,5A)£ 2 ^ 

Ja2(t) - A 



(9) 



to lowest order in A and SX. For a 1% strain, 5X/X = 0.01, 
and the same values for the other paramaters as in Eq. , 
we get a current density of ~ lOnA/mm. Eq. ([9]) does 
not contain A; this is because SX alone both breaks the 
rotation symmetry and tips the spins out of the a;y-plane. 



B. Physical process 

The appearance of the Berry curvature suggests a role 
of the anomalous velocity in generating the current. Such 
mechanisms have been discussed in the literature in the 
context of the CPGEfTil. [l6j . However, those mecha- 
nisms only work when the electric field changes slowly 



compared to the typical scattering time. The SSs of 
Bi2Se3 probably have lifetimes of tens of picoseconds, 
and thus, we are in the opposite limit when hui = O.leV, 
which corresponds to a time scale 10 3 times shorter. 

In this limit, the dc responses are a result of a preferen- 
tial absorption of the photon at one of the two momentum 
points for each pair of points (±p x ,p y ) related by to, as 
shown in Fig. QJl for p y = 0. According to the surface 
Hamiltonian ([7]), the spin vector S = ZH gets tipped out 
of the xy-plane for states that lie beyond the linear dis- 
persion regime, but the direction of the tipping is oppo- 
site for (p x ,p y ) and {—p x ,p y )- Thus, photons of helicity 
— 1, which can only raise (S z ) of an electron, are prefer- 
entially absorbed by the electrons that have (S z ) < in 
the ground state. The response, then, is determined by 
the properties of these electrons. Clearly, the process is 
helicity-dependent as reversing the helicity would cause 
electrons with (S z ) > to absorb the light preferentially. 

This is consistent with the requirement of a non-zero 
Berry curvature, which essentially amounts to the spin 
direction n having to be a three-dimensional vector. In 
the linear limit, where H = v F (p x a y — p y cr x ), the spin is 
entirely in-plane, and all the electrons absorb the incident 
light equally. 



C. Calculation in brief 

We now briefly outline the calculation of the helicity- 
dependent photocurrent. The detailed calculation can 
be found in Appendix [B] Readers only interested in our 
results may wish to skip this section. 

The Model: The Hamiltonian and relevant electric 
field (vector potential) perturbations for getting a direct 
current to second order in the electric field of the incident 
photon are 



H 
H' 

jo 

A x (t) + iA v (t) 



= \E p \n(p).a 



j x A x (t) 
dH 

dp a 



■jvMt) 



l e 



i{uj—ie)t 



(10) 

(ii) 

(12) 
(13) 



where A is the vector potential, z is assumed to be the 
surface normal, and e is a small positive number which 
ensures slow switch-on of the light. 

Quadratic response Theory: In general, the cur- 
rent along x to all orders in the perturbation H' is 



T* ( e'^'~ dt'H'(t') | ^ 



-ift^dtH't?) 



&>(*) 

(14) 

where T (T* ) denotes time-ordering (anti-time-ordering) 
and 0(t) = e Oe~ lHt . Terms first order in H' cannot 
give a direct current. The contribution to the current 
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from the second order terms can be written as 



symmetric tensor, the above traces can be written as 



t ti 

(j x )(t)= I dt' j dt"([\j x (t),H'(t')],H'(t")}) 



— oo — oo 
t U 



dt 1 / dt"x x ^(t,t',t")A a (t')Ap(t") (15) 



-oo — OO 



where a,P € {x,y}, Xx a p(t, t' , t") = Xxaff(0,t' - t,t" - 
t) = ([\jx,j a (t' - t)] ,jp(t" -t)}) = Xx«f}{t' - t,t" - t) 
due to time translational invariance, and the expectation 
value is over the ground state which has all states with 
E p < (>) filled (empty). For Hamiltonians of the form 
of Eq. (|10p , the expectation value of any traceless oper- 
ator O in the Fermi sea ground state can be written as a 
trace: 



Tr (HO) 



2 IK 



(16) 



This gives, 

/, . \ - Tr(H[[j x , j a (ti)} ,jp{t 2 )\) nr? s 
zapytiM) = - 2_j i 1 ') 



Eq. (H7| is the zero temperature limit of the finite 
temperature expression for the quadratic susceptibility 
proven in Ref. (l7j . 

Because of the mirror symmetry m, Xxap(ti, t 2 ) is non- 
vanishing only for a ^ j3. To get a direct current, we re- 
tain only the non-oscillating part of A x (t + ti)A y (t+tj) = 



A 0„2et 



sin (2u>t + Lo(ti + tj)) — sin (uj{ti — tj))]. Thus, 



dti j dt 2 | (Xxxy ~ Xxyx) (h,t 2 )x 

XJ 

(* 1+ * a )sin(w(t a -ti))) (18) 



-00 —00 



The Result: After carrying out the two time- 
integrals, we get the three currents mentioned in Eq. 
@. For clean samples at low temperatures, j a2 (t), 
which grows linearly with time, is expected to domi- 
nate. A general expression for this term is (in the units 
e = h = vp = 1 where vp is the Fermi velocity) 



ja 2 {t) 



iA§7rfsgn(oj) 
2J 2 



£*(M - 2\E p \)Tx(Hj x )Tr(H\j x ,j y }) 

(19) 



Using Eqs. (fTTJj) and ([P2")) and the Lie algebra of the 
Pauli matrices, [er^ix,-] = 2ieijk(^k where eijk is the anti- 



Tv(Hj x ) = 2\E p \v x (p) 
Tr(H[j x ,j y }) = A^n.f^-x^- 



(20) 



v dp x dp v t 
= U\E p \ 3 F(p) (21) 

Eqs. ([11]), (T2DD and flU) give our main result Eq. ©. 

IV. SPIN GENERATION 

Having understood the microscopic mechanism under- 
lying the generation of the photocurrent j a 2(t) , we won- 
der, next, whether such a population imbalance can lead 
to any other helicity-dependent macroscopic responses. 
Since each absorbed photon flips the 2-component of the 
spin of an electron, a net (S z ) is expected to be generated 
on the surface. 

The calculation of (S z ) is identical to that of jc pg e ■ 
The total (S z ) generated consists of the same three parts 
a s jcPGEi and the dominant part is 



%(*) = - 



7re 2 £g ftisgn(u;) 



^£(%|-2|K|H(p)F(p) 



(22) 

S z does not break the rotational symmetry of the surface, 
so we calculate S* 2 (t) directly for the threefold symmetric 
Hamiltonian ([7]) and obtain 



s z a2 (t) 



e 2 £ 2 {hwf\ 2 t 



A 



(23) 



For the same values of all the parameters as for j a2 (t) , we 
get S* 2 (t) ~ 10H, which means only ten electron spins are 
flipped over an area of ~ 1mm 2 . If, instead, we ignore 
the cubic corrections but assume magnetic ordering on 
the surface, so that H = vp(p x a y —p y cr x ) + Ma z , we get 



s z a2 (t) 



e 2 £ 2 M 2 t 
16(M 3 



(24) 



which again gives a rather small value of ~ 10h for 
M ~ 10K, a typical magnetic ordering transition tem- 
perature. However, the spin generated could be measur- 
able if one uses a pulsed laser of, say, Mega Watt power, 
and performs a time- resolved experiment. 



V. CONCLUSIONS 

In summary, we studied the CPGE on the surface of 
a TI at normal incidence, and applied the results to the 
(111) surface of Bi2Se3. If the rotational symmetry of 
the TI surface is broken by applying an in-plane mag- 
netic field or a strain, we predict an experimentally mea- 
surable direct photocurrent. A striking feature of this 
current is that it depends on the Berry curvature of the 
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electron bands. Such a dependence can be understood 
intuitively as a result of the incident photons getting ab- 
sorbed unequally by electrons of different momenta and 
hence, different average spins. The current grows linearly 
with time until a decay process equilibrates populations, 
which provides a way of determining the excited states 
lifetime. We also calculated the amount of dc helicity- 
dependent out-of-plane component of the electron spin 
generated. This does not require any rotational symme- 
try breaking; however, the numerical value is rather small 
with typical values of parameters. 

In the future, we hope to find a generalization of 
our results for oblique incidence. Experimentally, this 
is a very attractive way of breaking the rotational 
symmetry of the surface; indeed, such experiments have 
already been performed successfully on graphene[l9j. 
In graphene, helicity-dependent direct photocurrents 
have also been predicted by applying a dc bias(l8|. 
However, with a dc bias across a TI surface and ordinary 
continuous lasers, we find the current to be too low to 
be measurable. Finally, we also wonder whether the 
Berry curvature dependence of the helicity-dependent 
response to CP light survives for three- and higher-band 
models. If it does, it would be interesting to write such 
a model for semiconductor quantum wells such as GaAs 
and SiGe. It could also enable one to treat oblique 
incidence, by considering transitions to higher bands of 
different parities, because they are driven by the normal 
component of the electric field, E z . 
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tal inputs, and Ashvin Vishwanath and Yi Zhang for in- 
valuable feedback on the draft. 
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Appendix A: Proof of Berry curvature expression 

Here we show that the Berry curvature defined for 
Bloch electrons as 

F(p) = i ((d Pm u\dp y u) - (dp y u\d Pm u)) (Al) 

can be written as 

f(p) = n.(3 Pi nx5 Pj ii) (A2) 

for the band with energy | for Hamiltonians of the 
form H p = \Ep\n(p).(r. 

At momentum p, the Bloch state \u p ) with energy \E p \ 
is defined as the state whose spin is along n(p). Defining 
| f) as the state whose spin is along +z, \u p ) is obtained 
by performing the appropriate rotations, 

|ttp> = e -^(p) e ^ e (P) | t) (A3) 

where 8(p) and </>(p) are the polar angles that define 
n(p): 

n(p) = sin 9(p) cos <fi(p)x + sin 9(p) sin </>(p)y + cos 6(p)z 

(A4) 



Substituting Eq. (|A"3"]) in Eq. ([AT]) , one gets 

F(p) = sin0(p) (d p J(p)d p J(p) - d p J{p)d Py 9{p)) 

(A5) 

which, on using Eq. (|A4[) and some algebra, reduces to 
the required expression Eq. (|A2[) . 

Appendix B: current calculation for the cpge 

Here we explain the current-calculation of Sec. IHI Al 
in more detail and also state results for the parts of the 
current that we chose not to focus on there. 

As shown in Sec. MI CI the relevant susceptibility is 

X x ^(t,t\t") = -I^Tr^^W.r^')] ,/(*")]) 

= -E^^W^i)],/(* 2 )]) 

EE X Xaf) (t u t 2 ) (BI) 

where t\ = t' — t, t 2 = t" — t, and the non-vanishing 
components of x xa<3 are those for which a ^ j3. The 
non-oscillating part of the current, hence, is 

ti 

— oo — oo 

(X xxv (h,t 2 ) - X xyx (ti,t 2 )) e^ +t ^ sin (w(t 2 - t x )) 

(B2) 

Since jcPGEit) is an odd function of w. it reverses on 
reversing the polarization, as expected. 

The traces in the susceptibility expressions are calcu- 
lated by introducing a complete set of states in place of 
the identity several times. Thus, 

x xxv {hM) (B3) 

p nrnl 

^( El -E m)tl e -m-E n)t ^ Xnl X lm Y mn + C.C.j 

where X n i — (n\j x \m) etc. and the subscript p on E p 
has been dropped to enhance the readability. Similarly, 

X xyx (h,t 2 ) (B4) 
p p 

p nrnl 

(e^-^ XnlYlm _ e -m-^)t iYniXim ^ +c ^ 
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Substituting (|B3J) and (|B4J in (??), we get 

ti 

jcPGEit) = - J ^^z J dh I dt 2 e e ^ +t ^X (B5) 



— OO — OO 

J(E m -E n )t 2 



sin (ui(tt - t 2 j) ^2 sgn(E n )i 

p.nml 

{(e^-^* 1 - e -< E '- E ^) X nl X lm Y mn - 



where fRe stands for 'the real part of. Carrying out the 
the two time integrations gives 



A 2 e 2et 

jcPGE(t) = — Jm2222sgn(E n ): 

p nml 



(B6) 



1 



1 



E m - E n + u - ie E m - E n - to - ie 



Ei -En- 2ie 



E l -E m + 2ie 



where 3m stands for 'the imaginary part of. Using 

^(n 2 ^) = and *e(is=*) = h in the limit 

e 0, we get after some algebra, jcPGE(t) = j na +j a i + 



ja2(t), where (Tr denotes the trace) 

Al ^ io{io 2 - 12£j 



Jna 



E 



16 ^i\E p \ 3 {uj 2 -4^2)2 



Tr^gTr^bx,^]) (B7) 
comes from intraband processes and is constant in time, 



Jal 



E 



32 ^ 

p p 

Tr(H\j x ,\j x ,j y }}) 



(B8) 



is a result of an interband transition absorption as indi- 
cated by the <5-function in energy and is also constant in 
time, and 

, . Aiirt sgn(u;) v-^ n , ^ is 



Tr^^Tr^^]) 



El 



(B9) 



which also results from interband absorption and in- 
creases linearly in time. The last term was the main 
focus of our work. 
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